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ABSTRACT 


A model has been developed for the cooling of a quiescent hot layer by a 
sprinkler spray in a two-layer zone model (with no direct interaction of the 
spray with a fire plume), based on existing models of spray-induced flow and 
heat transfer to evaporating drops. In addition, existing models have been 
adapted to predict the penetration of the spray-induced flow below the layer 
interface and associated entrainment of lower-layer fluid into the upper 
layer. The cooling model is in good agreement with published results from 
sprinklered room-fire experiments, for which penetration of spray-induced flow 
into the lower layer and associated entrainment were calculated to be negli- 
gible. No published results have been found to test the validity of the 
adapted models of penetration and associated entrainment, which do not account 


for possible effects of the spray below the layer interface. 
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I 
INTRODUCTION 


The objective of this program has been to develop a submodel for sprin- 
kler/hot layer interactions to be incorporated in a two-layer zone model of 
compartment fires. 

In fire tests conducted within the large volume of the Factory Mutual 
Research Corporation test building in West Glocester, RI, the smoke layer 
(hot, upper layer) underneath the ceiling has been observed in numerous tests 
to remain well stratified during the active phase of a fire. However, smoke 
has been observed to penetrate into the clear air, especially during smoky 
fires. This penetration forms as pockets beneath sprinklers operating away 
from the immediate fire plume. It is clear that the smoke is pulled down by 
the action of the water droplets and returned to the smoke layer by buoyancy. 
A turbulent, churning motion is evident in the penetration pockets. Operation 
of the last sprinkler in a test usually signals that the sprinklers have 
gained control of the fire and gas temperatures under the ceiling are begin- 
ning to fall rapidly. Soon afterward smoke is transported completely to the 
floor by the operating sprinklers as buoyancy decays, and smoke is rapidly 
dispersed throughout the test volume. 

In the two-layer zone model of the objective, it is visualized that there 
is an upper layer of hot fire gases of uniform temperature under the ceiling 
of the compartment and a lower layer of relatively cool gases of uniform 
temperature extending to the floor. A single sprinkler operates near the 
ceiling within the upper layer. The spray entrains hot gases from the upper 
layer, producing a vertical jet toward the layer interface, Figure 1. The 
gases in the jet have increased water content and reduced temperature relative 
to the quiescent hot layer because of evaporation of water drops. After the 
jet penetrates the interface, it is surrounded by cool air and experiences an 
upward, buoyant force which eventually reverses the downward flow back into 
the upper layer, in the process entraining air from the lower layer, Figure 
2. The basic information required includes the evaporation rate of water and 
the mass transfer rate from the lower layer to the upper layer, from which the 
contribution to cooling of, and mass addition to the upper layer can be calcu- 


lated. 
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Figure 1. Entrainment of gas from hot layer by discharging sprinkler 
and formation of vertical jet. 
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Figure 2. Penetration of spray induced jet into lower layer, accompanied 
by entrainment from lower layer. 


FACTORY MUTUAL RESEARCH CORPORATION 


OTIN2.RU 


Entrainment of gas (air) into a water spray has been treated by Heskestad 
et al.(12), Penetration of a jet across a density interface has been inves- 
tigated by Turner‘3) , and entrainment of fluid of the penetrating jet into the 
layer from which it projects from the other side of the density interface, has 
been studied by Baines. (4) Finally, Yuen and Chen (9) have established heat- 
transfer correlations for evaporating liquid drops which make it possible to 
calculate evaporation and cooling rates. In this report the various pieces 
are adapted and assembled to produce the required prediction methods. 

For sprinklers operating near the fire source, one should ideally account 
for the direct interaction with the fire plume as well, including additional 
cooling and penetration of water drops to the base of the fire. The present 
work did not consider such effects. 

Elements of the problem under study have previously been investigated by 
others. Bullen ®) proposed a model for assessing the likelihood of a sprin- 
kler discharge causing a hot smoke layer to break up and form a plume down 
into the clear area, involving the downward drag of the aggregate of drops 
versus the buoyancy of the hot layer. Morgan! 7) proposed a tentative model 
for heat transfer to a sprinkler spray from a hot layer, based on heat trans- 
fer to a single water drop as it describes its trajectory, which he subse- 
quently compared with a few experimental results. (8) On the basis of Morgan's 
model, Heselden‘ 9) derived an empirical equation for the heat absorption rate 


of a sprinkler spray, also quoted by Hinkley 10) | 
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jm 
ELEMENTS OF MODEL 


2.1 ENTRAINMENT INTO WATER SPRAY 

A one-dimensional theory for entrainment of air into water sprays 
discharging vertically downward in a quiescent, isothermal environment has 
been derived and experimentally verified by Heskestad et al. (12), The 
essential features of the theory include a prescribed spray boundary and 
momentum exchange between the water drops and the entrained air. 

The theory leads to the following set of coupled differential equations 


for the drop velocity and air velocity in the spray. 


4-1 -1 3/2 
gae/ox* = (€/6) u* - —.6u* u* - u*) 1) 
D £/8) uh hs ( 
du®/dx® = - u®x*7! , Bx" > (us ux)7! (us = ux)3/4 (2) 
with the initial conditions that at x* = (2 tane)~ |, otahar =a 
where: 
us = aba ek nondimensional drop velocity 
Se — U/UD9) nondimensional air velocity 
¥F cs x/D,,, nondimensional axial distance 
Uno = initial drop velocity 
D, = effective initial spray diameter 
The initial drop velocity is calculated from: 
Ti PROUL Wis qineya oy (3) 
po M W 
where Cy, is a momentum coefficient, defined as the initial momentum of the 
spray to the momentum of the water flow in the spray nozzle; v is the 
volumetric discharge rate of water from the nozzle; and D is the nozzle 
diameter. Furthermore, the effective initial spray diameter is defined as: 
y -1/2 
ae get ee (4) 
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The parameters y, 6 and &€ in Eqs. (1) and (2) are nondimensional, 


defined according to: 


Y= 5 oe (tan 5) * Me 
8 =..(3/4, Byaseuppimn Cle in e/* ay 7e - 
E =463/4) Bv'/"p9 ~' gapicanis 73/2 pi/4 7 


This is a convenient set for independent study of effects of spray angle, 6 
(via y), nozzle discharge velocity, u, (via 8), and nozzle diameter, D (via 
&). Other symbols include p,,, the water density; p, the air (gas) density; and 
g, the acceleration of gravity. The symbol r is a drop size parameter in the 


correlation equation: 
d = r(D/u,)°/3 (8) 


A value r = 0.084 mes~°/3 is typical of standard sprinklers when d is inter- 
preted as the volume-mean diameter in the spray. The value Cy = 0.4 is repre- 
sentative of momentum coefficients of standard sprinklers‘), B is a constant 
of proportionality related to aerodynamic drag of drops, given as 12.6, (2) 
Table I presents tabulated values of nondimensional entrainment flow for 
a spray angle (at room temperature) of © = 120° (y = 0.03445+103), taken from 
Reference 2. (Results are also presented for © = 90°, 60°, 30° in Ref. 2.) 
Nondimensional entrainment flow, Q*, is presented as a function of nondimen- 
sional vertical distance from the spray device, x*, at various combinations of 


8 and €, where: 
O&e= MEAVE (9) 


Here, Ve is the total volumetric gas flow rate over a cross section of the 
water spray. 

Table II presents nondimensional particle (drop) velocities associated 
with the results listed in Table I. These calculations have not been reported 
previously; they have been transcribed from the same computer printout from 
which Table I was constructed. (2) 
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TABLE I 


NUMERICAL SOLUTIONS FOR ENTRAINMENT IN CONICAL sprays?) 


-3 


Q* at x* indicated 


voto *pei0s eet io, eo) HOS 80160320 KD. 1280. 4.2560 
103445 1 7.5 239 712 2,060 5,710 14,600 32,400 72,400 176,000 465,000 
(120°) 15 239 713 2,070 5,770 15,100 37,900 99,300 272,000 761,000 
30 243 735 2,180 6,400 20,300 53,100 151,000 431,009 1,230,000 
60 283 899 2,820 8,730 26,700 80,000 236,000 681,000 1,950,000 
2 7.5 308 919 2,620 6,860 15,400 31,600 64,600 136,000 312,000 
15 308 919 2,620 6,870 15,500 33,200 75,400 189,000 516,000 
30 308 922 2,640 7,000 17,000 42,600 114,000 320,000 912,000 
60 314% 956 2,830 8,150 23,100 65,600 187,000 535,000 1,520,000 
4 7.5 384 1,140 3,120 7,390 15,500 31,500 63,400 128,000 260,000 
15 384. «1,140 3,120 7,390 15,500 31,600 64,700 137,000 317,000 
30 384 1,140 3,120 7,410 15,700 33,600 77,500 200,000 562,000 
60 385 1,140 3,150 7,690 18,100 46,400 128,000 367,000 1,070,000 
8 7.5 462 1,340 3,410 7,500 15,500 31,500 63,400 127,000 255,000 
15 462 1,340 3,410 7,500 15,500 31,500 63,400 128,000 260,000 
30 462 1,340 3,410 7,500 15,500 31,600 64,800 138,000 321,000 
60 462 1,340 3,410 7,530 15,800 33,900 79,100 199,000 599,000 
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TABLE II 
NONDIMENSIONAL DROP VELOCITIES ASSOCIATED 
WITH THE CALCULATIONS IN TABLE I 


u* at x*®* indicated 


yo10-2 8-10" F-10" 10 420 fo 80 160 320 640  72c0mmNeneD 

03445 1. 7.5 .943 °° 874 .735. 491.177. 9050308 <.0366» . O306NNOREATS 

(120°) 15.048 "883 .755.05..531 922668) MeN TAt eenlS7oc Geen 
30. 1.019 . 1.022), 1.016, -°.493)° 4.055) i920. 8e-902 dene BOn meee 
bo TB 231 5.27 5.53 5.54 5.5K 


sy 


_ 
Ls) 
Wo) 
oO 
wW 
foe) 
oO 
A 
(oa) 
=, 


2 7.5 .905 .788 .569 .260 .0685 .0295 .0155 .0091 .0061 
15 .906 .789 .570 .263 .0758 .0415 .0291 .0232 .0201 
30 910 .798 .589 .304 .149 121 .110 . 103 .098 
60 .981 .934 .843 .719 .636 .603 .586 $525 .567 


4 7.5 .852. .673 . .389 «141. .0574' “20272 9960134); | (O0CRI ae 
Ne 652 .673 .389 .142 .0578 .0280 .0146 .0083 .005H 
30 652 .674 .391 .145 .0636 =—.0364 .0244 8 .0187 .0157 
60 O01 a. 6.0035. SUOme aC 1Ee eth 0926 .0805 .0730 .0679 


8 7.5 ..186. .550,. »2714, 1166 0544 990265 632013127 -OGGl Meee 
if) 786. .§50,. «271s. o116— «20544 950266 e520132j;¢ .OOGTRgoee 
30 786, ~550.— »2Tlay «116% sOS4% ge0272 eae0141 .0079  .0049 
60 786 .551 272 .118 ~~ 0589" "20333". 0215™ {Ot GOMiaammas 
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The results in Tables I and II, presented as functions of x* for a given 
spray angle © (120°), can alternatively be presented as functions of the non- 
dimensional radius, R* = R/D,, where R = x tan (0/2). If the tabulations for 
other spray angles are consulted?) | it is found that Q* and u* as functions 
of R* are not very sensitive to the spray angle. Hence, Tables I and II can 
be used to predict entrainment and drop velocity for any spray angle likely to 
be encountered with automatic sprinklers, tying predictions to spray radius, 
rather than spray angle. 

The particular association of spray angle © with the parameter y in 
Tables I and II has been calculated from Eq. (5), assuming room temperature 
airacp 221.20 kg/m?) . At other temperatures the relation between y and © is 
obviously different, as is the relation between the spray radius R and the 
spray angle 6 discussed in the preceding paragraph. If room temperature 
conditions are identified with subscript "o", it may be verified using Eq. (5) 
that the equation for spray radius in the preceding paragraph is generalized 
to elevated gas temperatures in the spray, Ts, according to: 


R= x (TT) 
Sa40 


tan (6/2) (10) 

The isothermal theory has to be assessed for effects of buoyancy forces 
in the application to sprinkler sprays discharging into a hot layer. Assuming 
a certain temperature of the gas phase in the spray, lower than the surrounding 
hot gas, one may calculate ratios of buoyancy forces to drag forces of the 
water drops, both tending to propel the entrained gas, building on elements of 
the existing theory. A specific case, considered previously in its isothermal 
version!) , has been calculated. This case involved a water discharge rate 
(Vo) of 1.89°1073 m3/s (30 gpm), an initial drop velocity (Uno) of 15 m/s, and 
a volume mean drop diameter (d) of 0.4 mm. At the position in the conical 
spray where the cross sectional area (A) was 0.66 n°, the drop velocity (u,) 
was 9 m/s and the air velocity (u,) was 3.7 m/s. As a perturbation, we assume 
a surrounding hot gas (upper layer) temperature of ay = 200°C and a gas tem- 
perature in the spray of T, = 100°C. For this particular case, the ratio of 
buoyancy forces to drag forces is estimated at 0.10, sufficiently small so as 
not to be dominant in the process. Effects of other conditions may be 


assessed based on the following proportionality: 
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Buoyancy forces A ie 0.85 see “ 3/2 
Drag forces ee (7) a hina oe (11) 
V Ss u p s 


From this proportionality it is clear that buoyancy forces can become impor- 
tant in a continuously expanding spray. A tenfold increase in the cross- 
sectional area of the spray, not unrealistic in sprinkler applications, can 
make the buoyancy forces comparable to the drag forces. Effects on gas and 
drop velocities in the spray can be assessed from Eqs. (1) and (2), together 
with calculated results in Tables I and II and similar calculations for other 
values of y. It is found for the case of a 120° spray (sprinkler-like) at x* 
= 80, using typical values of 8 (2+107°) and & (30-1077), that including 
buoyancy forces equal to the drag forces increases the local entrainment flow 
rate by 28 percent and the local drop velocity by 18 percent. Effects of 
buoyancy forces are neglected in this investigation but might be included in 
future..studies. 

Other assumptions which are necessary in order to apply the existing 
entrainment predictions are 1) that evaporation of water is not sufficient to 
have a signifjgcant effect on the number density and size of the droplets 
propelling the entrained flow, and 2) that gas motion in the upper layer is 
not sufficient to invalidate the assumption that the surroundings of the spray 
are quiescent. With respect to the former assumption, it is sufficient to 
show (a posteriori) that the mass evaporation rate of water in the upper layer 
is small compared to the total discharge rate of water. With respect to the 
latter assumption, it is observed that significant motion in the upper layer 
will be confined to the ceiling jet generated by the fire plume, where velo- 
cities will typically be an order of magnitude smaller than the local drop 
velocities, which are expected to be close to their discharge value. 

It is convenient to express the values of the spray parameters in Eqs. 
(1) - (3) in terms of their values at room temperature (20°C). If this is 
done it is found: 

5 


-0.1 
8, (T,/T,) (12) 


hoe) 
i" 


-0.15 
AS (13) 


m 


10 
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where By and bo are calculated using room temperature. Instead of expressing 
y in terms of Tap it is more useful to relate associated values tan (6/2) and 
tan (60/2): 


tan(6/2) = (Tate tan (6/2) (14) 


Hence, the value 6 = 120° in Tables I and II refers to somewhat larger angles 
at elevated temperatures, as already implied in Eq. (10). 


2.2 PENETRATION OF SPRAY INDUCED FLOW INTO LOWER LAYER 

The penetration depth of a jet against opposing buoyancy forces into a 
fluid of different density has been studied by Turner. (3) In application to 
the problem at hand, the jet can be considered to be the flow entrained from 
the upper layer by the water spray, directed downward, whose characteristics 
are described by the jet velocity, width and excess temperature relative to 
the lower layer. 

Turner studied the upward rise (penetration) of salt-water jets injected 
into fresh water. Using dimensional arguments, he predicted that the upward 
rise would be a linear function of a parameter formed from the momentum flux 
and buoyancy flux of the jet, assuming a point source. He plotted his experi- 
mental values of rise against this parameter for three nozzle sizes and 
demonstrated good linearity, as predicted. However, for small values of rise, 
the experimental points appear to deviate slightly from the adopted straight 
line fit, as if there were effects of finite nozzle sizes. 

We decided to replot Turner's data on the assumption that effects of 
finite nozzle size could not be ruled out. The relationship sought was the 
nondimensional jet rise or penetration, PAS versus the nozzle Froude number, 


Fr, where 6p is the jet rise, R is the nozzle radius, and: 

Fr = u/({g(40/0,)R] (15) 
In Eq. (15), u is the average nozzle velocity; Ap is the density difference, 
salt water versus fresh water; and Py is the salt water density. Figure 3 


presents the results of nondimensional rise versus Froude number for the three 


nozzle radii of the experiments, based on the data in Turner's Figure 5(3) 


11 
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R=7.0: mm 
R=4.8 mm 
R= 35.5 ,.mm 
REGRESSION 


Fr 


Figure 3. Replot of Turner’s data on upward rise (penetration) of salt- 
water jets injected into fresh water, using nondimensional coordinates. 
Fr, defined in the text, is the Froude number at the injection nozzle. 


12 
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(and corrections of the abscissa of that figure, communicated to us by 
Turner‘11)), It is apparent, in this version of the data, that a minimum 
Froude number on the order of unity is required to produce jet penetration. 
The curve is a regression fit of the form y = a (x-b)”, specifically: 


0.870 (16) 


srg = 3.69 (Fr - 1.1) 
In adaptation to the interaction of sprinklers with hot layers, Eq. (16) 


may be written: 

6 /R, = 3.69 (Fr - 1,1)°°87 (17) 
where R; is the spray radius at the level of the interface and 

Fr = u,,/(g(sT,/T,) Relies (18) 


Here, u,; is the (uniform) gas velocity calculated in the spray at the level 


of the setae 

Although this adaptation of Turner's work (3) is appealing, it must be 
considered tentative. In the adaptation, the jet flow in the hot layer is 
generated by a spray whose effects have been ignored below the layer inter- 
face. The spray below the interface will entrain surrounding gas, which may 
affect the jet penetration. Furthermore, the spray produces a downward drag 
on the penetration, which is opposed by an upward buoyancy force on the 
penetrated fluid. Both of these forces are very nearly proportional to the 
penetration depth and, therefore, it is not clear what the net effect on the 


penetration depth will be. 


2.3 ENTRAINMENT OF LOWER LAYER FLUID INTO UPPER LAYER 

Baines(4) has investigated turbulent entrainment at a density interface 
like that between the upper and lower layer in the present problem. This 
entrainment refers to the entrainment of lower layer gas into the intruding 
jet from the upper layer, which reverses back into the upper layer. 

In Baines' case, the density interface was produced in a tank by placing 


a layer of fresh water over a layer of salt water. Jet experiments were per- 
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formed with the aid of a small pipe producing a jet of fresh water against the 
density interface (from above), and plume experiments were performed substi- 
tuting salt water for fresh water in the jet discharge. Dye was introduced 
into the jet to observe the entrainment process. Quoting from Baines‘4); 


"As the plume impinged on the interface it was observed to entrain 
fluid through its end over a relatively small area. This entrained 
fluid mixed with the plume fluid within a small volume, i.e., about 
the volume of a large eddy. Upon impact with the interface the 
large rotating lumps of fluid, i.e., the largest eddies, stopped 


their forward motion and moved out along the interface." 


Baines found that a nondimensional volumetric entrainment rate could be 
expressed as a function of a Froude number, where the characteristic length, 
velocity and buoyancy were all associated with the plume or jet as it impinged 
on the density interface. His entrainment results can be expressed: 


V 
= errs (19) 


= 5.13-10 G 


2 
ted 
where ve is the volumetric entrainment rate from the heavier layer; w, is the 
local centerline velocity (at impingement) of the plume or jet; b, is the 


local radius to 0.368 times the centerline velocity in the plume or jet, 
assuming a Gaussian profile; and the Froude number is defined: 


ws 


Ena 


72 (20) 
(@) 


In the expression for the Froude number, p54 is the density of the heavier 
fluid (salt water); p, is the local density on the plume or jet centerline; 
and p,, is the density of the lighter fluid* (fresh water). 


Although Baines' Froude number is defined with Pao, it is believed that 
p, is a more appropriate density at the location of p,, in Eq. (20). 
Baines' correlation, Eq. (19), would not be affected significantly since 
Pao and 0, were of similar magnitudes in the experiments. 
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Baines' Froude number, Eq. (20) is based on variables of a Gaussian pro- 
file, while the Froude number in Eq. (15), used to correlate Turner's data in 
Figure 1, is based on plug flow in a tube. The two Froude numbers can be 
related, requiring that the two flow profiles have the same momentum flux, 


buoyancy flux, and total mass flow rate, so that; 


us 4 Wi bp/o, = ‘ (p5- ,)/0.,3 R = ut b (21) 
with the consequence: 
Fre = 23/4 Ep (22) 
Baines(4) noted the depth of penetration into the heavy fluid layer, 


Stating it “was barely perceptible at small Frig] but for Frig] > 2, the plume 
or jet penetrated deeply." Using Eq. (22), the inequality Frag > 2 becomes Fr 
> 1.2, and Baines' observation then appears quite consistent with the inter- 
pretation of Turner's data‘3) in Figure 3. 

In terms of Fr, Eq. (19) becomes: 


Ve 3 
<Sape 0.244 Fr 5 (23) 
uR 


This form is compatible with the form of the induced flow predictions dis- 
cussed in Section 2.1, which produce average values for the spray cross 
section. However, one more refinement is needed in application to the fire 
problem in recognition of the possibility that a large difference in density 
may exist between the upper and lower layers. In analogy with the work of 


(12) 


Ricou and Spalding on local entrainment rates in jets of different density 


than the surroundings, the mass entrainment rate across the interface is 


expected to obey: 
: 1/2 
nm, * R (Mo, ) (24) 


where M is the local momentum flux in the impinging jet. Substituting for M, 


one obtains: 
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m, = R-u(o.e,) (25) 
Hence, a Suitable nondimensional entrainment rate is defined as: 
n, : m /0, f Ne 
rn VEPeF Tieise 
2 1/2 2,°s 2-8 
uR’(p.p,) uR C ) uR* (==) 
L s 
Consequently, Eq. (23) is generalized to: 
v Tomi/2 
oS exoruh (22m Meer (26) 
2 iT 
ear s 
ih EB 


where Me = M/P—s the volumetric entrainment rate from the lower layer, and u 
and R have been replaced by us; and R; for application to the spray/hot layer 
problem in which case Fr is interpreted according to Eq. (18). 

As cautioned at the end of the preceding section, there may possibly be 
effects of the spray traversing the penetrated fluid and of air entrained into 
the spray below the layer interface, which have not been accounted for in this 


treatment. 


2.4 HEAT TRANSFER TO EVAPORATING LIQUID DROPS 

Yuen and Chen‘?) measured heat transfer to simulated (porous sphere) 
water and methanol droplets in a vertical hot air tunnel and established a 
relationship for the Nusselt number in terms of the Reynolds number, Prandtl 
number, and mass transfer number. The following expression can be derived 
from their work for the convective heat transfer coefficient for evaporating 
water drops, using the approximation that the dynamic viscosity of air is 
proportional to the 0.7 power of temperature: 


1/2 1/3 -0.05 
patti b. p(T/T,) (u_- u,) d C Ho Tp 
pe (TAT oa k aS 
h Ousluec fe) (27) 
3 1 -0.85 


1.19 
Cleat (T, - T.) ) 
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where: 
h = convective heat transfer coefficient for a drop 
a= density of air at standard temperature and pressure (STP) 
Ty = room temperature (20°C) 
T, = air temperature 
uy = particle (drop) velocity 
u, = air velocity 
d = drop diameter 
lee A dynamic viscosity of air at room temperature 
Tp = lie + Tp)/2; film temperature 
T, = temperature of upper layer 
Tp = drop temperature (equal to wet bulb temperature) 
Cy = specific heat of air at constant pressure 
K, = thermal conductivity of air at room temperature 
r eb 100. K-19 a constant 
‘i = 273K, a reference temperature 


The film temperatures, Tp, require knowledge of the wet bulb tempera- 


tures, Tp» given by Yuen and Chen‘ 13) as a function of the air temperature, 


T,: 


Temi sc) 
100 
200 
400 
600 
800 

1000 


Th (°C) 


ce 
50 
70 
85 
90 
20) 


These values can be used to establish the following expression for the film 


temperature: 


Tp(K) = 0.543 (TAK) + 245) 
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This model for the convective heat transfer coefficient assumes that 
water drops contributing significantly to evaporation and cooling in the hot 
layer are in approximate thermodynamic equilibrium with the environment, i.e., 


all heat transferred to the drops is absorbed in vaporizing the drops. 


2.5 HEAT BALANCE IN UPPER LAYER 
Equation (27) becomes, after substitution for Tr from Eq. (28) and 
various constants, as well as substituting the subscript "s" for "a" to refer 


specifically to the gas phase in the spray: 


T-2,| Leal 0.40 
ye cecil bul es fe) ; 1/2 
5.16-10 ~+ 4.55-10 sy roa | (u,- u,)d 


0.85 (29) 
| ate 


T 
-4 1.19 _O 
(1 + 1.16-10 (T.- 273) ) E 


h( kW/m°-K) = 


The heat transfer per drop in the spray is: 
2 
had Gre ~ 7 
The drop generation rate is: 
. 3 _ . 3 
Vi /( ad /6) = 6 V/(nd ) 
The drop transit time in the upper layer is: 
(6 - 6.)/u, 


Consequently, the heat absorbed by the water spray can be calculated from: 


V (6 - 6_) 
0) cane (har alee eee 
s Siete i xs q3 Uy 
V Lees 6.) 
Q. = 6 a nn cates (sete 7) h (30) 


where h is given by Eq. (29). 
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The heat removed from upper layer gas entrained in the spray is equal to 
the sum of heat absorbed by the water drops and the heat absorbed by air 


entrained from the lower layer: 


" w Sues : 2 
an Coal Tee Tu ER i eae ot h + D, Cy (T. T,) (31) 


s and the drop diameter d. Equation 


(31) is solved for the temperature in the gas phase in the spray, Tg, assuming 


with h given by Eq. (29) in terms of Up» u 


that this temperature is constant throughout the spray. The upper layer ten- 
perature is taken as the value in the absence of the spray. The quantities 
Up» Us, V, and d are calculated with the aid of the methods in Section 2.1, 
and the mass flow entrained from the lower layer, mo) is calculated as des- 
cribed in Section 2.3. The heat absorbed by the spray is then obtained from 
Eq. (30) and the mass addition rate of water vapor to the upper layer is 
calculated from: 

nate Q./L (32) 


where L,, is the latent heat of vaporization of water (2260 kJ/kg). 


The drop and gas velocities in the spray, u, and us, vary with the 


distance from the sprinkler. In the heat Aa calculations, values 
averaged over the distance form the sprinkler to the layer interface are 
employed. 

It has been assumed so far that the drops in the spray are uniform in 
size, or at least can be represented by their volume mean diameter. Whereas 
the volume mean diameter can be justified for entrainment calculations('»2) | 
we have not been able to show that the volume mean diameter is the appropriate 
effective diameter for cooling calculations. However, the effective diameter 


for cooling, d., might well be proportional to the volume mean diameter, d: 


c 


defzird (33) 


where ~ is the constant of proportionality, as yet not known. 
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In many cases the spray induced flow in the upper layer will not pene- 
trate significantly into the lower layer, according to the model discussed in 
Section 2.2. However, in other cases significant penetration may occur. In 
cases with significant penetration depth, one may assume that the heat trans- 
fer coefficient, gas temperature and relative velocity between drops and gas 
are the same as their respective averages in the spray above the interface. 
The average drop transit time through the penetration pocket is assumed 
as n 6_/u_, where n is a constant of perhaps 0.5, 6p is the (maximum) depth of 
the penetration, and Up is the average drop velocity above the interface. 


Then the revised form of Eq. (31) would have the following substitution for 
foo 6): (6 - b5 + 7 6). 
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1G Uy! 
VERIFICATION AND CALIBRATION OF MODEL 


3.1 AVAILABLE EXPERIMENTS 

A series of 25 sprinklered room fires have been reported by You et al. 
(14,15) which are available for testing the cooling predictions of the 
model. The tests involved three (steady) heat release rates, three sizes of 
geometrically similar sprinklers, and three sprinkler discharge rates. 

A plan view of the test room! 14) is reproduced in Figure 4a, showing a 
fire source located near one end of the test room and a floor-to-ceiling 
opening at the opposite end. The (pendent) test sprinkler was mounted under 
the ceiling, near the end of the room with the fire source, with its deflector 
0.10 m from the ceiling. The three available geometrically similar sprinklers 
had nozzle diameters of 11.1, 8.36 and 6.94 mm. A side view of the fire 


source | 14) 


is reproduced in Figure 4b, a heptane spray fire generated by a set 
of four oil-burner nozzles. The nozzle assembly was shielded from direct 
sprinkler water impingement with a 3.1 mm thick steel plate. The total heat 
release rate of a fire and the convective heat loss through the room opening 
were determined by calorimetry of the fire gases rising from the room open- 
ing. With the additional measurements of heat losses to ceilings and walls 
and radiative losses through the room opening, it was possible to determine 
the heat absorption rate of the sprinkler spray. 

The 25 tests conducted by You et ai. (14,15) are listed in Table III; the 
full table has not been reported previously but has generously been made 
available by the authors '6) for the purpose of this study. The table lists 
the sprinkler diameter(D); the discharge rate (Vo); time from ignition when 
the measurements began; the actual heat release rate (Q,)5 the convective heat 
loss rate (Q,)5 heat losses to ceiling and walls, and radiative losses to 
floor and room opening (Q,)5 the heat absorption rate of the water spray (Q,); 
the distance from the ceiling to the neutral plane in the room opening (h,); 
the centerplane bulk gas temperature rise of the hot flow in the door 
opening! '5) (OT); and the ambient temperature (T_). 

Three freeburn tests were conducted, i.e., Test 13 at a heat release rate 
of 134 kW, Test 17 at 246 kW, and Test 21 at 396 kW. It may be verified that 


the bulk temperature rise, aT scales nearly as Gee for these tests, as 
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(a) Plan view of test roon. 
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(b) Heptane spray fire source. 


Figure 4. Test facility employed by You et al sata 
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TABLE II] 
conpiTions‘'®) cor tests conpucTED BY you ET aL. (14,15) 


; Beginning 
D Me ttt %3 Q& Q, Qs, hy, aT, 
Test No. (mm) (Liter/min) Period (s) (kW) (kW) kW) (Kw) (m) (2) 
1 6.94 45.4 300 150 24.0 22.9 103niaee 0770 30 
2 6.94 68.2 300 148 HeOT gn 10.7 124.7 0.76 1 
3 6.94 45.4 300 314 65.0 41.1 207 29° “21.00 =| 
4 6.94 68.9 300 30601152256) 823076 252.6 0.91 25 
5 6.94 45.4 300 4o4 = 122.0 += 63.4 308.6 1.07 87 
6 6.94 68.2 300 486 8 =58.0 43.0 385.0 0.96 45 
| 8.36 40.1 300 149 53.6 27.0 68.4 0.92 49 
§ 6.36 68.2 300 140 15.1 20.3 104.6 0.80 24 
9 &.36 45.4 300 304 101.1 46.5 156.4.) 91.05 1 
10 8.36 68.2 340 316 57.4 38.4 220.2 0.94 50 
11 8.36 45.4 303 44y7 = 162.6 += 80.9 20355912 F. O77 P1417 
12 8.36 68.2 300 472) 104.3, S44 213834 [1205 68 
13 11.10 0 250 134° «91.0 844.3 0 0.87 85 
14 ch heat] us 4 300 135 1320 esi s.t 30.9 0.98 68 
15 eee 68 300 1377 ©40.0 421.1 15; Soe Vente a1 
16 etd 98.4 360 136 16.0 14.1 105.9 0.91 23 
17 17.10 0 280 246 173.0 90.6 0 1/05 > 136 
18 11.10 45 360 262) 164.0 60.3 57.1 1201e8 118 
19 11.10 68 300 282 103.0 44.2 134. 1.03 77 
20 11.10 98.4 360 285 54.0 28.1 202.9 1.08 52 
21 11.10 0 140 396 228.0 145.6 0 tee) Aller 
22 Ti10 us .h 300 418 220.0 106.6 91.4" 1.16" 161 
23 11.10 68.2 300 460 158.0 72.4 229-0 wala lo, tic 
24 11.10 68.2 300 474 157.0 76.4 240.6 1.24 113 
25 11.10 98.4 300 453. 96.0 =46.2 310.8 1.22 69 


= timed from ignition. 
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(17) (16) 


expected based on known scaling rules. Vertical temperature profiles 
in the centerplane of the room opening for these tests are shown in Figure 5, 
presented as AT/AT, versus 2/h,, the latter being the distance from the ceil- 
ing normalized with the distance of the neutral plane from the ceiling. A 

fair approximation to a constant-temperature upper layer is evident, as 
represented by the dashed lines which define approximately the same integrated 
area as the actual temperature profile. 

Among the sprinklered tests listed in Table III, a number of them produced 
temperature profiles in the room opening which can be considered to be of 
approximately constant temperature, as presented in Figure 6616), These tests 
included all the tests for which Q,/Q, was 0.55 or less. A number of tests 


(16) produced approximately linearly 


With a ratio Q./Q, of 0.65 or greater 
decreasing AT/AT, with increasing 2/h,, reaching zero typically near 2/h,=1.5. 

It was decided to limit analysis of sprinklered tests for the most part 
to cases with approximate two-layer temperature profiles, i.e., Tests 7, 9, 11 
for the 8.36 mm nozzles and Tests 14, 15, 16, 18, 19, 20, 22, 23, 25 for the 
11.10 mm nozzle. Among these, Tests 16 and 20 (with Q./Q, values of 0.77 and 
0.71) apparently did not produce approximate two-layer profiles. Test 24 was 
not included since it was a repeat of Test 23, with very similar results. 


of their three geo- 
(18) 


You et al. measured the volume median diameters,* d., 


metrically similar sprinklers at different water pressures and found them 


to be correlated by an equation like Eq. (4): 
: 2/3 | 
d= 0.136 (D/u,) (34) 
Assuming a ratio of volume-median to volume-mean diameters of 2.1, as 
previously reported for a Rockwood T-4 nozzle(®) | the coefficient r in Eq. 


(4), pertaining to the volume-mean diameter of the nozzles used in the 


experiments, becomes: 
r = 0.065 mes~¢/3 (35) 
The volume median diameter divides the drop population into equal volume 
fractions of smaller and larger drops. 
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Figure 5. Freeburn temperature profiles in centerplane of room opening ‘1° 


in experiments reported by You et al. 
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Figure 6. Temperature profiles in centerplane of room opening for tests with 
it Ae = 0.55 and smaller in experiments reported by You et gia Aase? 
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This value is slightly smaller than the value 0.084 mes~°/3 quoted after Eq. 
(4) as typical of standard sprinklers. 


3.2 VERIFICATION OF COOLING MODEL 


Verification of the cooling model will be investigated by comparing model 
predictions against the experiments summarized in the preceding section, in 
the following manner. For each experiment, the effective diameter for cooling 


introduced in conjunction with Eq. (33), d., will be calculated from the model 


, 
so that the predicted cooling rate sel, cts measured cooling rate. It must 
then be shown that the associated fraction ¢ = d./d [Eq. (33)] is approximately 
constant among the various experiments, where the volume-mean diameter, d, is 
calculated from Eq. (8). Using the value for r in Eq. (35), a realistic cali- 
bration for ¢ is also established in this procedure. 

First the properties of the hot layer in a two-layer representation is 
established. It is observed in Table III that the depths of the neutral 

plane, h,, of the tests selected in the preceding section range from 1.0 to 

1.2 m, including the freeburn tests. From Figures 5 and 6, the depth of the 
upper layer in a two-layer representation may be taken as z/h, = 1, corres- 
ponding to actual depths, by , in the range 1.0-1.2 m. An average depth gy 

= 1.1m, will be assumed for all the selected tests. Referring to Figures 5 
and 6 again, it is seen that a representative upper layer temperature corres- 


ponds to AT//AT, = 0.9; actual values of temperature rise in the upper layer, 


st are then iene as 0.9 AT), where AT, is listed in Table III. For 
each sprinklered test, a value of the freeburn temperature at the heat release 
rate of the sprinklered test is needed. This temperature is determined using 
the freeburn test closest in heat-release rate to the sprinklered test, then 


He, as justified in Section 3.1. 


adjusting the value aT» assuming aT, « Q. 
Calculation of the spray induced flow and drop velocities in the upper 
layer requires information on r, B and Cy (Section 2.1). The nominal values 
0.084 mes~°/3 | 12.6 and 0.4 are selected, respectively. Although the nominal 
value for r, 0.084, is slightly larger than the best value available for the 


experiments* in Eq. (35), the difference is of little consequence for the 


The result in Eq. (35) was not available at the time the calculations were 
performed. 
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entrained flow calculations‘! , The typical value Cy = 0.4 is assumed because 
an experimental value is not available for the tested sprinklers. Other 
information required includes 1) the nozzle discharge velocity, u,, calculated 
from the discharge rate and nozzle cross-sectional area; 2) the effective 

w? calculated from Eq. (7); and 3) the initial drop 
velocity, u,., calculated from Eq. (9). 


initial spray diameter, D 


To proceed further, Test 23 is selected to illustrate procedure. 
Referring to Table III, this test involved a nozzle diameter D = 0.0111 m and 
a sprinkler discharge rate Ve = 68.2 £/min = 1.136-1079 m/s. The following 


quantities are calculated: 


0 /(aD*/4) = 1.136-1073/ (x 0.0111°/4) = 11.7 m/s 


Cc 
iT 


RAal/20 


-1/2 
wW M 4 


- 0.0111 = 0.0176 m 


oO 
t 


0. 


. 2 
use = Cy [V /(xD /4)) = 0.4-11.7 = 4.68 m/s 
The parameter g may first be evaluated at room temperature, using Eq. (2) and 
B= 12.6, v = 1.490107? m/s, p = 1.2 kg/m, p, = 1000 kg/m>, Cy = 0.4, r = 
0.084, u, = 11.7 m/s: 


1.54°107°, 


ww 
iT) 


The parameter € at room temperature follows from Eq. (3): 
E> = 45.8610" 


At the gas temperature in the spray, 1 which is larger than room tempera- 
ture, the values of 8 and € are smaller than B. and €, according to Eqs. tie) 
and (13). It is assumed for the present purpose that Ts can be set equal to 
the upper layer temperature, Tih evaluated as described earlier in this sec- 
tion: T, = 126°C = 399K. Hence (T,/T,)~°-'> = 0.954, so that: 

= 1.47°107¢ 

= 43.7°1074 


~~ wo 
i] i) 
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Table IV lists Q* and uy" as functions of x* for the values of 6 and € 
just calculated, interpolated from Tables I and II and pertaining to a room 
temperature spray angle of 120°. The nondimensional parameters are converted 
to physical variables using Eqs. (5), (6) and (8). The spray radius corres- 
ponding to a given axial location x in the spray depends on the gas tempera- 
ture in the spray, Tes which is not known. However, as a first trial value 
the upper layer temperature determined for Test 23, T = 399K, is selected. 

The spray radius can then be calculated from Eq. (10) (using ©,/2 = 60°) and 
the gas velocity in the spray is determined from: 
u_ = V_/aR? (36) 
s Ss 

It appears from spray distribution measurements( 16) that the spray radius 
of the test sprinklers at the depth of the layer interface below the sprinklers, 
1.1 - 0.1 = 1.0m, is well represented by the radius of a 120° spray cone, 
which is 1.73 m. Hence, a location x is sought in the spray, represented in 
Table IV, where the spray radius is 1.73 m at the trial value of Ts. According 
to Eq. (10), this radius occurs at x = 0.856 m. The values of Up and u, for 
this location are determined by interpolation between x = 0.707 m and 1.408 m, 
and the associated value of ve is obtained from uc(aR). 

Penetration of the spray-induced jet into the lower layer and entrainment 
of lower-layer fluid into the upper layer may now be investigated, using the 
relations in Sections 2.2 and 2.3. There is no penetration and the entrain- 
ment rate of lower layer fluid is three orders of magnitude smaller than the 


flow rate, Ve , associated with the spray induced jet at the layer interface. 


i 
With insignificant entrainment from the lower layer, the heat balance in 


Eq. (31) can be written in the convenient form: 


Bilaate! sa CN tle (37) 
where b = Vogt o/Ty) CS (38) 
a= 6 (V/d) ((6,-8.)/u, Jh (39) 


29 


FACTORY MUTUAL RESEARCH CORPORATION 


OTIN2.RU 


TABLE IV 
TRIAL CALCULATION OF GAS AND DROP VELOCITIES 
IN SPRINKLER SPRAY OF TEST 23 
(g = 1.47°1072; & = 43.7°1074) 


x* Q* Te x(m) Ns (=H u (5) Trial T,(K) R(m) Us ) 
10 284 1.05 OC176 «0.322 4 95 399 07805 0.813 

20 869 1.02 0.352 0.987 4.81 399 Oras 0.621 

40 2550 0.99 O77 07mnee Oo 4.67 399 1.42 0.456 

mit. ~ = 0.856 4.03 4 42 399 lays 0.429 

80 7340 0.75 1.408 8.34 3.54 399 2.84 0.329 

ay ety Sh - a ee aapces «39 

s2 u 0 (T/T )C 4 .03-1.2 (293/474)-1.00 


a One Oma U aan 


(T,5 is essentially equal to first trial value of T,, 399K; hence, T, = 397 K.) 
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From Eq. (30) it is observed: 
a(T,-T)) = Q, (40) 
Combining Eqs. (37) and (40): 
b(T, -T.)= Q. (41) 
from which: 
fet, 8.72 (42) 


Eq. (42), with substitution for b from Eq. (38), has been used at the bottom 
of Table IV to test the trial value of i (in terms of Tyo). Values of iby and 
Q, measured in the experiment are employed, together with the value of ver 
resulting from the trial value of Ts; 4.03 m3/s. It turns out that Ts2 is 
essentially equal to the trial value in this case, so T, = To5 = 397K is the 
final calculated value. If T,5 had not matched the trial value, it would have 
been employed as a second trial value and the procedure continued until 
convergence occurred. 

The drop size which produces the correct cooling, Qs: according to the 


model can now be investigated. From Eq. (40): 


de: (43) 


With substitution from Eq. (42), one obtains: 


Q. 


L 2 2a perc) S 
ae - Q./b 


Q 


eet (44) 


(T,-1,) ~ Q/(V, (067, /T) C,) 


For Test 23, Q. = 230 kW; T, = 474K; upper layer temperature (experimental 
substitute for T,) = 126°C and the associated Tp = 39°C [see list following 
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Eq. (27)]; and Voi = 4,03 m°/s. The value a = 2.64 kW/K is obtained, which is 
substituted in Eq. (39), together with the expression for h from Eq. (29) to 
determine the effective drop diameter for cooling, d, [see Eq. (33)]. Values 
of Up and u, are taken as averages of values ljsted in Table IV over the 
distance from the sprinkler to the pseudodistance, x = 0.856 m (true 

distance eur = 1.00 m), where the calculation spray attains a radius of 
1.73 m. Calculation of averages requires values at the location of the sprin- 
kler (x = x* = 0); consistent with the theory(2) , the initial values of both 
Uy, and u, are set equal to Uno: 
4.70 m/s and 1.02 m/s, respectively.) With the appropriate substitutions in 


(Average values calculated for Up and u, are 


Eq. (29), the effective drop size for cooling is calculated (iteratively) as 
d, = 0.336°1073 m for Test 23. This drop size may be compared to the volume 
mean diameter for the conditions of Test 23, calculated from Eq. (4) using the 
value of reingEq. (35)',\i.ea ied = 0.627°1073 m. Using the definition in Eq. 


(33)5 
ee d/d = 0.54 


The other selected tests have been treated similarly. As in the case of 
Test 23, no penetration of the spray-induced jet into the lower layer was ever 
evident, and the calculated entrainment rate of lower layer fluid into the 
upper layer was always less than one percent of Vays i.e., the calculated flow 
rate associated with the spray induced jet in the upper layer at the layer 
interface. 

Figure 7 presents all the ratios ¢ = d,/d calculated, plotted as a func- 
tion of the heat release rate. The various data points are identified by 
their associated test numbers. The plot covers essentially three heat release 
rates, three discharge rates of the sprinklers, and two sprinkler sizes. It 
is observed that there is a slight decline in ¢ with heat release rate, con- 
trary to ideal behavior and the reason for which is not clear. The arithmetic 


mean of all determinations is: 
c = 0.70. (45) 


This level is indicated in the figure. 


32 


FACTORY MUTUAL RESEARCH CORPORATION 


OT1N2.RU 


O 100 200 3500 400 500 
Q ( kW ) 


Figure 7. Ratios [ = d/d indicated by experiments of You et ag AR 


plotted as a function of the heat release rate. 
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The calculations represented in Figure 7 assume, it is recalled, that the 
drop temperature is the wet bulb temperature corresponding to the calculated 
gas temperature in the spray. Calculations assuming a drop temperature equal 
to the wet bulb temperature at the freeburn temperature have shown a slightly 
greater decrease of ¢ with Q, whereas calculations assuming a drop temperature 
equal to the ambient temperature have indicated a dependence of ¢ on Q similar 
to, Figure: 7. 

Assuming ¢ constant at the value in Eq. (45), one may now calculate the 
associated sprinkler cooling rates in the selected experiments and compare 
these with the directly measured rates. The effective drop diameter for 
cooling is calculated from Eqs. (45), (33), (8) and (35). Next, the value of 
a in Eq. (39) is calculated with the aid of Eq. (29), assuming T, is unchanged 
from the previous calculations. From Eq. (37): 


bT - aT 
ON, onsen ak 58) (46) 


Ss a+b 
where b is unchanged from the previous calculations. Equation (46) is used to 
check the original assumption for Ts. If the new and original values differ, 
the new value for T; is used in a new calculation of a. The procedure is 
continued until the value for T, from Eq. (46) matches the assumption to 
calculate a. Finally, the cooling rate is calculated from Eq. (40), using the 
established values of a and Ty. The results of these calculations are 
presented in Figure 8 as Calculated Q, versus Measured Q.. 

The procedure described above is acceptable since the final Tg does not 
differ greatly from the original assumption for Ty. The two values may differ 
Significantly in a normal calculation of cooling rate in the intended applica- 
tion of this model. In that case, a trial value for vs is first assumed, such 
as the freeburn temperature in the upper layer, associated values g and £& 
evaluated, and then the equivalent of a table like Table IV prepared, yielding 
a first value for Vues This value for vies allows b and a to be calculated and 
the first trial value of Ts to be checked. If the new and the trial values 
differ (as they probably always will), the new value is substituted as a 
second trial value and the procedure repeated a sufficient number of times 


that Eq. (46) confirms the last trial value. 
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Figure 8. Comparison of calculated and experimental Se spray cooling 


rates. (Test numbers indicated.) 
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IV 
PENETRATION AND ENTRAINMENT BEHAVIOR 


In Sections 2.2 and 2.3 it was concluded that applications of Turner's 
work (3) to the penetration behavior of the spray-induced jet into the lower 
layer and of Baines' work 4) to the associated entrainment of lower layer 
fluid into the upper layer had to be considered tentative. The reason for 
this caution is that there may possibly be significant effects of the spray 
traversing the penetrated fluid and of air entrained into the spray below the 
layer interface, which have not been accounted for in the model. Nevertheless, 
a few calculations with the predictive relations developed are instructive. 

In Section 3.2, based on applications of Turner's(3) and Baines! (4) 
results, none of the experiments employed to verify the cooling model produced 
jet penetrations or significant entrainment from the lower layer. In this 
section, penetration and entrainment behavior is investigated using various 
upper layer depths and gas temperatures within the spray. 

It is assumed that a 1/2 in. standard sprinkler is operating at the 
ceiling level (D = 0.437 in. = 11.1 mm), discharging at Vv, = 1.89-107> m/s 
(30 gpm) and producing a 120° spray over the range of ceiling layer depths 
considered. Other assumptions include: Cy = O Mies r= 502005 m/s?’ 3, E = Oe); 
T, = 20°C. Associated quantities include D, = 0.0176 m; u, = 19.6 m/s; Uno 
7.84 m/s; B, = 1.99°107¢: and &) = 4.57°1073, Entrainment tables were esta- 


blished, similar to Table IV, using appropriate values of 8 and & for each 
value of T, selected. Equation (17) was used to calculate penetration depth 
and Eq. (26) was used to calculate entrainment rates from the lower layer. 
The results of the penetration calculations are presented in Figure 9 as 
penetration depth Bs versus the difference in temperature between the gas in 
the spray and the lower layer, T, - T,. Various smoke layer depths are 
represented: 0.35 m, 0.71 m, and 1.41 m. Except for the smallest depth, 
temperature differences are rather small for significant penetrations. Once 
penetration begins to occur as temperature differences decrease, the pene- 
tration depth increases rapidly with further decreases. (A penetration depth 
extending to the floor might be interpreted as a breakdown of the two layer 
model, with smoke dispersal throughout the enclosure volume.) Figure 10 


presents associated volumetric entrainment rates from the lower layer 
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Figure 9. Influence of the temperature difference between gas in spray and 
lower layer on jet penetration depth beneath layer interface for various 
layer depths, Ou (1/2 in. standard sprinkler discharging 1.89 2/s at 
ceiling level.) 
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Volumetric entrainment rates from lower layer associated 


with Figure 9. 
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(Vo), including a brief table for reference showing entrained flow rates in 
the spray at the interface level (V.). Significant entrainment rates from the 
lower layer develop as the temperature differences between the gas in the 


spray and the lower layer is reduced. 
No data have been found to test the predictions in Figures 9 and 10. 
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V 
MODEL PARAMETERS 


The following input parameters are indicated: 


D sprinkler nozzle diameter (m) 
volumetric sprinkler discharge rate (m3 s7!) 
R; spray radius at interface between upper and lower layers (m) 


é depth of upper layer (m) 


: distance of sprinkler below ceiling (m) 

T, upper layer temperature in absence of spray (K) 

1% lower layer temperature (K) 

Cu Momentum coefficient of sprinkler, with a typical value of 0.4. 

r volume-mean drop size parameter, with a typical value of 0.084 mes~°/3 
C ratio of effective-cooling drop diameter to volume-mean drop 


diameter in spray, with a typical value of 0.70 


Fixed input parameters include: 

v kinematic viscosity of air at room temperature (1.49°107> me s7') 
p density of air at STP (1.20 kg m3) 

p density of liquid water (1000 kg m7 3) 
g acceleration of gravity (9.81 m s~*) 
T 


reference temperature (273 K) 


The following output parameters have been recognized: 


Ve volumetric spray-induced flow rate of gas in upper layer (m3 s7') 

us spray induced downward gas velocity in upper layer (m s7!) 

Uy drop velocity in upper layer (m s7!) 

T, average gas temperature within spray in upper layer (K) 

Q, heat absorption rate of water spray (kW) 

m, mass entrainment rate of gas from lower layer into upper _jayer, generated 
by spray-induced flow impinging on layer interface (kg s'') 

Dey mass evaporation rate of water (kg s7!) 

Ps penetration Hepes of spray induced gas flow below layer interface (m) 
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VI 
CONCLUSIONS 


A model has been developed for the cooling of a fire-heated hot layer in 
a room by a sprinkler spray, based on existing models of spray generated 
flow and convective heat transfer coefficients for evaporating drops. 

The model assumes upper (hot) and lower (cool) layers of uniform tempera- 


tures and no direct spray interaction with a fire plume. 


The cooling model has been augmented with models for penetration of the 
spray-induced hot-layer jet into the lower layer and associated entrain- 
ment of lower-layer fluid into the upper layer, based on previous studies 
with momentum driven jets and buoyant plumes. It is cautioned that these 
models are tentative because of possible effects of the spray on the 
penetrating fluid. According to these models, significant penetration 
and entrainment from the lower layer occurs only at low hot layer tem- 
peratures (compared to typical temperature ratings of sprinklers), with 


greatest effects for shallow hot layers. 


The cooling model has been tested against published results from 
sprinklered room-fire experiments, where the hot layer depth was about 
1.0 m below the sprinkler. Good agreement was observed under various 
conditions, which included two sprinkler sizes, three water discharge 
rates and three heat release rates, although a segment of the sprinkler 
spray intercepted the fire plume. The effective cooling diameter of the 


drops in the spray was indicated at 70% of the volume-mean diameter. 


Calculated penetration depths and associated entrainment rates from the 
lower layer were insignificant in all experiments referred to in the 
preceding conclusion. No data were available from the experiments to 
verify these calculations. Furthermore, no other data have been found to 
test the penetration and associated entrainment predictions in 


sprinkler/hot layer situations. 


yy 
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SYMBOLS 

A Cross sectional area of spray 
a defined in Eq. (39) 
B constant in Eqs. (2) and (3) 
b defined in Eq. (38) 
by local radius to 0.368 times centerline velocity in plume or jeC 
Cy ratio of initial spray momentum to flow momentum in nozzle 
Cy specific heat of air at constant pressure 

nozzle diameter 
he effective initial spray diameter; see Eq. (7) 

drop diameter, or volume mean drop diameter in polydisperse spray 
Fr Froude number; see Eq. (15) 
Frg Froude number, see Eq. (20) 
f 1.16°1074K7 1-19 a constant 
g acceleration of gravity 
h convective heat transfer coefficient 
ta distance from ceiling to neutral plane 
Ky thermal conductivity of air at room temperature 
Li latent heat of vaporization of water 
M local momentum flux in impinging jet 
m, mass entrainment rate from lower layer into upper layer 
mi, mass addition rate of water vapor to upper layer 
Q, actual heat release rate 
Q, convective heat release rate 
Q, heat absorption rate of spray 
OX WE nondimensional entrained flow in spray 
R spray or nozzle radius 
R* R/D,,, nondimensional spray radius 
R; spray radius at interface level 

drop size parameter defined in Eq. (8) 
Te air temperature 


Tr (Tit T)/e, film temperature 
Te lower layer temperature 


ak room temperature (20°C) 
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drop temperature 

273K, a reference temperature 
gas temperature within spray 
temperature of upper layer 
ambient temperature 
temperature rise 

bulk temperature rise of fire gases leaving room 
velocity 

air velocity 

nozzle discharge velocity 
drop velocity 

initial drop velocity 

gas velocity within spray 

u. at interface level 


Ss 


Un /Un9) nondimensional drop velocity 


volumetric entrainment rate from lower (heavier) layer into upper 


(lighter) layer 


volumetric gas flow rate in spray 

Ve at interface level 

volumetric water discharge rate of spray nozzle 
centerline velocity of plume or jet (in Baines' model) 
distance below spray nozzle 

x/D,,, nondimensional distance below spray nozzle 
distance below ceiling 

viscous interaction parameter; see Eq. (2) 

8 evaluated at room temperature 

geometric parameter, see Eq. (1) 

y evaluated at room temperature 

penetration depth of jet 

distance of sprinkler below ceiling 

depth of upper layer 

constant of proportionality, see Eq. (33) 
constant 

conical spray angle 


dynamic viscosity of air at room temperature 
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kinematic viscosity of air 


gravity 


parameter; see Eq. (3) 


—& evaluated at room temperature 


density 
density 
density 
density 
density 
density 
density 
density 
density 


of gas (essentially air) 
of air at STP 

of jet fluid 

of heavier fluid 

of lighter fluid 

of lower layer 

of gas within spray 

of water 


difference 
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